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Abstract 



We consider a quantum mechanical system of N bosons with relativistic dispersion interact- 
ing through a mean field Coulomb potential (attractive or repulsive). We choose the initial wave 
function to describe a condensate, where the N bosons are all in the same one-particle state. 
Starting from the iV-body Schrodinger equation, we prove that, in the limit N — * oc, the time 
evolution of the one-particle density is governed by the relativistic nonlinear Hartree equation. 
This equation is used to describe the dynamics of boson stars (Chandrasekhar theory). The 
corresponding static problem was rigorously solved in |10j . 

o 
in 

1 Introduction 

We consider a system of N identical bosons with relativistic dispersion relation and with a mean 
field Coulomb interaction. The system is described on the Hilbert space l? s (M. 3N , dx) , the subspace 
of L 2 (M 3Ar ,dx) containing all functions symmetric with respect to permutations. The Hamiltonian 
of the system is given by 

N , N 

Here we use units with H = m = 1 (where m denotes the mass of the bosons). The Hamiltonian Hjy 
defines a mean field interaction among the bosons because the coupling constant is proportional to 
1 /N: with this scaling the kinetic and potential part of the energy are typically of the same order. 
This condition is necessary in order to have a mean field description of the system in the limit of 
large N. 

The constant A can be positive or negative, corresponding to repulsive and attractive interaction. 
In the case of repulsive interaction we will have no restriction on the value of A. The choice of a 
negative A, which corresponds to an attractive Coulomb potential, leads to a Chandrasekhar theory 
of boson stars, where general relativity effects are neglected. In this case we need to impose the 
condition A > — 4/ir. To understand why this assumption is necessary, note that, if A > — 4/7T, the 
Hamiltonian Hjy, with domain D(Hn) = H 1 (M. SN ), is self adjoint and stable, in the sense that the 
ground state energy of Hn divided by the number of particle N, is bounded below, uniformly in 
N. On the other hand, for A < — 4/V, the Hamiltonian Hn is unstable: the ground state energy 
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per particle diverges to — oo as N — > oo. This phenomenon is usually referred to as the collapse of 
the system: the energy is minimized by letting the particles closer and closer, because the increase 
of the kinetic energy (due to the localization of the particles) is not enough, if A < — A/tt, to 
compensate for the decrease of the potential energy. These results were derived in ^U], where the 
authors prove that the energy per particle of the ground state of the system is given, in the limit 
N — > oo, by the minimization of a one-particle energy functional: 

lim — = inf £(ip,~0) 

N^oo N ¥>:|MI=1 

(1.2) 



inf ( fdx |(1 - A) 1 ^ 4 ip(x)\ 2 + £ fdx(V 
3: llvll =1 \J 2 J 



Here we defined V(x) = The r.h.s. is negative infinity if A < — 4/V. The existence of a 

critical coupling A cr ;t = — 4/7r is due to the fact that the kinetic energy, which behaves, for large 
momenta, as |V|, and the potential scale in the same way. 

In this paper we are interested in the dynamics generated by Hn in the limit N — > oo. From 
(|1.2j) one can expect that the macroscopic dynamics of the system are described, in the limit 
N — > oo, by the one particle nonlinear relativistic Hartree equation 

id m (x) = = (1 - &) 1/2 Vt(x) + X(V * \<f>t\ 2 )(x)<pt(x) . (1.3) 

dip 

Some important properties of this equation, such as the global well-posedness for A > A cr i t , are 
proven in [SJ. 

To formulate the convergence towards the nonlinear Hartree equation 1)1. 3J) more precisely, we 
define next the marginal distributions of an iV-particle wave function ^jy, and we investigate their 
time evolution. Given an iV-particle wave function V(x) € L 2 S (R 3N ,dx) with ||^|| = 1, we define 
the corresponding density matrix 77V = |V'A r )(V ; -/v I as the orthogonal projection onto ip. The kernel 
of 7jv is 

7at(x;x') = VaKx^a^x') . 
More generally, a density matrix 7^ is a non-negative trace class operator on L^(lR 3Ar ) with Tr 77V = 



(k) 

1. For k = 1, . . . N, the A;-particle marginal distribution of 77V, denoted by "y N , is defined by taking 

(k) 

the partial trace over the last N — k variables of 77^, that is, the kernel of "y N is given by 



Here and henceforth we use the notation 

x = (x 1 , . . . ,x N ), x k = (xi, . . . ,x k ), x N _ k = (x k+1 ,x k+2 , ■ ■ ■ ,x N ) 

with Xj S M 3 , for all j = 1,...N, and analogously for the primed variables. Since Tr77\r = 1, it 
follows immediately that Tr 7^ = 1 for all k = 1, . . . , N. 

The time evolution of the system is governed by the Schrodinger equation 

id t i>N,t = H N ip N)t (1.4) 
or, equivalently, by the Heisenberg equation 

idtjN,t = [HN,lN,t] 
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for the dynamics of 77V • From the Schrodinger equation we can also derive a hierarchy of N equa- 
tions, usually called the BBGKY hierarchy, describing the evolution of the marginal distributions 

(k) 

j\ri for k = 1, . . . , N. Using the permutation symmetry we find, for k = 1, . . . , N, 

k 



^7g = E[( 1 -^) 1/2 ^S] + ^ E [v( 



\Xi - X 



3=1 



l<i<j<k 



;i.b) 



(fc+1) 

t 



3=1 



(Recall that V(x) = 1 ). Here we use the convention 7!^ = 0, if k = N + 1. Moreover we 

denote by Tr^+i the partial trace over the variable Xfc+i- In terms of the kernels 7^1 (x^x^) the 
last equation can be written as 



ft7$(x*;xi) = E ((! - A ^ /2 " ( X " A,0 1/2 )7^(x fe ;x' fc ) 



i/2^(*), 



3=1 



l<i<j<fc 



7 ^ P 

+ A(l - — ) E / dx fe+ i - x fc+ i) - F(a^ - x fc+ i)) 7^| 1) (x fe , x fc+1 ;Xfc,a? fc+ i) . 

3=1 

In the limit JV — > 00 the BBGKY hierarchy ()1.5|) formally converges to the infinite hierarchy of 
equations 



•H (fe) - E [(!- A -J "■-/ 



3=1 



3=1 



(fc+i) 



(1.6) 



for all k > 1. It is easy to check that, if the initial data is factorized, that is if 

k 

7£ ) (x fc ;x / fc ) = Y{<p(x j )Tp(x' j ) 



3=1 

for all k > 1, then the infinite hierarchy (|1.6ft has the solution 

k 

7f ) (x fc ;x' fc ) =Y[(f>t(xj)vi(x , j ) 
3=1 

where tpt(%) is the solution of the non-linear one-particle Hartree equation Q1.3|) with initial data 
<£>t=0 = V 9 (f° r A > A CT it = — 4/7r, Ifl.MJ) is known to have a unique global solution in the space 
H m / 2 (R 3 ), for every m > 1, see (§])■ Hence, if we consider a factorized initial wave function 
■0Ar(x) = ri^Li l f{ x i)i if we fix A; > 1 and f € R, and if we denote by 7^ the /c-particle marginal 
distribution associated with the time evolution il>N,t of tpN, then we can expect that, in a suitable 
weak topology, 



TN,t ( x fc ; 4 ) -» Tt ( x fc ; x fc) = II & to ) ^ (4 ) for ^ ~* 00 

3=1 



1.7) 
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where cpt is the solution of the Hartree equation Q1.3|) . The aim of this paper is to give a rigorous 
proof of this statement. 

The first result of this type was proven in [7] for non-relativistic dispersions and for smooth 
potentials. This work was generalized to bounded potentials in |13| . In |S] the authors show the 
convergence (|1.7(l . for a non-relativistic dispersion, and for any integrable potential: they use the 
formalism of second quantization and they need the initial state to be coherent (and thus the number 
of particles cannot be fixed). In [I] the convergence (|1.7|) was proven for bosons with non-relativistic 
dispersion interacting through a Coulomb potential. Partial results for the non-relativistic Coulomb 
case were also established in PQ. In j2], a joint work with L. Erdos and H.-T. Yau, we consider a 
system of N non-relativistic bosons, interacting, in the mean-field scaling, through an iV-dependent 
potential V N (x) = N^V{N p x), with < (3 < 3/5, which converges to a delta-function in the limit 
N — > oo: for this potential we prove the convergence of solutions of the finite hierarchy ()1.5|) (with 
V(x) replaced by V^{x)) to solutions of the infinite hierarchy (jl.6|) (with V replaced by S(x)). A 
recent overview of rigorous results and open problems concerning the nonlinear Hartree equation 
and its derivation as the mean field limit of large bosonic systems can be found in 5 . 

Note that the non-relativistic case considered in all these works is technically easier than the 
relativistic case we are considering in the present paper: the presence of a quadratic dispersion 
(the non-relativistic kinetic energy of a boson is given by the negative Laplacian, quadratic in the 
momentum) makes the control of the Coulomb singularity much simpler than in the relativistic case, 
where the kinetic energy is only linear in the momentum (for large momenta, (1 + p 2 ) 1//2 ~ \p\). 
Hence, although our general strategy is the same as in we need here more refined estimates. 

Next we explain the methods we use to prove (|1.7|) in some details. The main technical tool in 
our analysis is an a-priori estimate (Theorem 14. 4|) . which guarantees a certain smoothness of the 
iV-body wave function (and thus of the corresponding marginal densities), uniformly in N and t. 
As in J3j, we derive the a-priori estimate from energy estimates, which control powers of the energy 
in terms of the corresponding powers of the kinetic energy. In order to prove these energy estimates 
we need to introduce an iV-dependent cutoff in the Hamiltonian. For e > we define the cutoff 
Hamiltonian 

N N 1 

H N = V(l - Aj) 1 / 2 + A Y , } — r . (1.8) 

^ J ^ \xi- xA+ eN' 1 K ' 

j=l Kj 

In Hn the Coulomb singularity has been regularized on the length scale |x| ~ eN^ 1 . In Propo- 
sition we then prove the following operator bound (valid on the space L 2 (M 3Ar ) of functions 
invariant w.r.t. permutations) for powers of the Hamiltonian H^: 

H% > C k N k (1 - A 1 ) 1 / 2 (l - A 2 ) 1/2 ... (1 - Afc) 1 / 2 (1.9) 

for all k > 1, and for all iV large enough (depending on k and on the cutoff e). Hence, if the initial 
state ip N £ L 2 s (R m ) is such that 

(iP N ,H k N ip N )<C k N k (1.10) 

for all k > 0, then, using the conservation of the energy, from (|1,9[) we can derive bounds on higher 
derivatives of the solution i[) e N t of the Schrodinger equation 

idt^t = HN^N,t (1-H) 

with regularized interaction, and with initial data ip% t=0 = ipN- Unfortunately, a typical factorized 
wave function ^iv(x) = Ylf^^^j) does not satisfies (|1.1U() (unless ip G C^°(IR 3 ), an assumption 
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we want to avoid). Therefore, the introduction of an additional cutoff in the initial one-particle 
wave function (p is required. Given k > and a one particle wave function ip € -ff 1 (3R 3 ) , we define 
its regularized version by 

ip K = exp(— k\p\/N)<p with p = — iV . 

Then the regularized iV-particles wave function t/>^(x) = 11^=1 f K { x j) can be proven to satisfy 
(jl.lUj) for every k > 1 and for all N large enough (depending on k and k). This follows from 
Proposition where we derive an upper bound for the expectation of powers of the Hamiltonian. 
Hence if we denote by ip% K t the evolution of the regularized initial data ip 1 ^ generated by Hn, then 
(jl.9j) implies that 

(1055, (1 - A fc )V^«) < C k (1.12) 

for all k > 1 and for all iV large enough (depending on k, e and n). For fixed e, k > we denote 
by Tpft = {7^*}fcLi the family of marginal distributions associated with the wave function V'jvt- 

Using the bound (|1.12|) we can prove the compactness of the sequence T^t with respect to an 
appropriate weak topology. Moreover we can show that any limit point = {7^\}fc>i of T^t 

satisfies the infinite hierarchy ()1.6|) with initial value 7^\ = o( x fc; x fc) = Tlj=i Pfaj^Wj) (the cutoffs 
e and k, which regularize the interaction and the initial data, disappear in the limit N — > oo) and 
that the solution of (|1,6|) is unique (in an appropriate space). This implies that 

k 

for every k > 1, where <p>% is the solution of the non-linear Hartree equation ()1.3|) with initial data 
(ft=o = V 9 - Finally, we prove that the difference between the physical evolution (generated by the 
Hamiltonian H^) of the initial wave function t/>tv( x ) = Ylj=i<p( x j) an d the modified evolution 
(generated by Hjv) of the approximate initial wave function ip^ converges to zero, for e and k 
converging to zero, uniformly in N (for N large enough). This concludes the proof of (|1.7[) . 

The paper is organized as follows. In Section |^1 we introduce some Banach spaces of density 
matrices, and we equip them with weak topologies, useful to take the limit N — > oo. In Section|31we 
state our main result and give a sketch of its proof. In Section 0] we prove the energy estimates for 
the regularized Hamiltonian Hn, and we use them to derive the a-priori bound (|1.12j) . In Section 
we show the compactness of the sequence of marginals 7^ associated with the wave function ip e ^ K t 
with respect to the weak topology. In SectionElwe demonstrate that any limit point of the marginal 
densities Fjyj = {^\}k>i satisfies the infinite hierarchy (|1.6|) . In Section[3we show the uniqueness 
of the solution of (ITT)]) . In Section El we remove the cutoffs, letting e,k — > 0. Finally, in Section 
we collect some important technical lemmas, used throughout the paper. 

Notations. Throughout the paper, we will use the notation Sj = (1 +p|) 1//4 = (1 — Aj) 1//4 for 
integer j > 1. Moreover, with an abuse of notation, we will denote by \\tp\\ the L 2 -norm of the 
function ip and by ||^4|| the operator norm of the the operator A. We will use the symbol Tr^ to 
denote the partial trace over the fc'th variable. In general we will use the symbol Tr to denote the 
trace over all variables involved: sometimes, instead of Tr, we will use the symbol Tr^ fc ^ to stress the 
total number of variables over which the trace is taken. In general, we will denote by C a universal 
constant, depending, possibly, only on the coupling constant A and on the i? 1 -norm of the initial 
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one-particle wave function ip. If C also depends on other quantities, we will usually stress them 
explicitly. 

Acknowledgements. We are grateful to H.-T. Yau and L. Erdos for suggesting this problem to 
us, and for very useful discussions. We thank G. M. Graf for helpful suggestions. B. S. would also 
like to thank J. Frdhlich, S. Schwarz and E. Lenzmann for stimulating discussions related to the 
relativistic Hartree equation. 

2 Banach Spaces of Density Matrices 

For k > 1, we denote by C\ and by K, k the space of trace class and, respectively, of compact 
operators on the fc-particle Hilbert space L 2 (M 3fc ,dx fc ). It is a well known fact (see, for example, 
Theorem VI. 26 in that (£L || • ||i) = (fC k , || . ||)*, where || . ||i denotes the trace norm, and || . || 
the operator norm. This induces a weak* topology on C\. We define C 1 as the direct sum over 
k > 1 of the spaces C\, that is 

^ = 04 = {r = {7 ( %>i : 7 (fc) e 4, V k > 1} . 

k>l 

We equip the space C l with the product of the weak* topologies on C\ for all k > 1. A sequence 

Tat = {7^ }fc>i converges to T = {7^}fc>i in C 1 for N — > oo if and only if jffi converges to 7^ 
with respect to the weak* topology of C\ for all k > 1. 

To prove our main theorem (Theorem I3.1JI . we also need to define a different topology, which 

(k) 

makes use of the smoothness of the marginal densities 7]^ and allows us to pass the smoothness 
to the limit N — > 00. For 7W G C\, we define the norm 

||7 (fc) || Wfc =Ti\S 1 ...S kl ^S k ...S 1 \ 
and the corresponding Banach space 

n k = { 7 ^ eel-. \\^\\ Hk < 00}. 

Recall that we use the notation Sj = (1 +p|) 1//4 = (1 — A^) 1 / 4 . Moreover we define the space of 
operators 

A k = {T^ =Sl.. S k K^S k ...St: e /C fc } 

with the norm 

\\T {h) \U k = ll^r 1 • • • s^T^s^ 1 . . . sr 1 !! 

where | . | denotes the operator norm. Then we have, for every k > 1, 

CHk, II • Httfc) = (A, II • IU fc )* (2-1) 

and thus we have a weak* topology on the space H k . The proof of (|2.1|) is identical to the proof 
of Lemma 3.1 in [1]. Note that the weak* topology on TL k is stronger than the weak* topology on 
the larger space £\: if a sequence 7]^ converges to with respect to the weak* topology of TC k , 
then it also converges with respect to the weak* topology of C\. 

We will prove that the densities 7^ associated with the wave function ip% K t (defined by the 
regularized Schrodinger equation (jl.lljl ). are continuous functions of t with values in TL k (continuous 
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w.r.t. the weak* topology of Hk)- In order to prove the continuity of limits of 7^1, we want to 
invoke the Arzela-Ascoli Theorem, which is only applicable to metric spaces. Fortunately, since 
the space Ak is separable (which follows from the fact that the space K-k of compact operators is 
separable), it turns out that the weak* topology of Hk, when restricted to the unit ball of Hk, 
is metrizable. In fact, because of the separability of Ak, we can fix a dense countable subset of 
the unit ball of Ak- we denote it by {Ji^}i>i £ Ak, with ||<//^||.a,, = 1 for all i > 1. Using the 
operators J^' we define the following metric on TCk'. for >y( k \^( k ) g Hk we set 

oo 

P,(7 (fc) ,7 (fc) ) == E 2_i | Tr ^ (7 W "7 {fc) )| • (2-2) 

Then the topology induced by the metric pk{-,-) and the weak* topology are equivalent on the 
unit ball of Hk (see ^2]) Theorem 3.16). In other words, a uniformly bounded sequence 7^ £ TCk 
converges to E TCk with respect to the weak* topology, if and only if pk(^\l^) —> for 
N -> 00. 

For a fixed T > we consider the space C([0,T],Hk) of functions of t € [0, T] with values 
in Hk, which are continuous with respect to the metric pk (for uniformly bounded functions, the 
continuity with respect to pk is equivalent to continuity with respect to the weak* topology of Hk)- 
On the space C([0, T], Hk) we define the metric pk by 

p fc (7 (fc) (-),7 (fe) (-))= sup p fc (7 (fe) (i),7 (fe) (*))- (2-3) 
te[o,T] 

Finally, we define the space H as the direct sum over k > 1 of the spaces Hk , that is 

H = ($H k = {T = {7 (fc) } fc >i : 7 (fc) G Wfc, VA; > 1} , 
fc>i 

and, for a fixed T > 0, we consider the space 

C([O,r],H) = 0C([O,T],H l ), 
fc>l 

equipped with the product of the topologies induced by the metric jo^ on C([0, T], Hk)- That is, 

for T^t = {7yv,l}fc>i and r * = {7t fc) }A:>i in C([0, T],H), we have -> T t for TV -> 00 if and only 
if, for every fixed A; > 1, 

pkiiNtiTt ) = sup Pk{r Nt ,r t ) -* 

te[o,T] 

for iV — > 00. 



3 Main Result 

We are now ready to formulate our main theorem. We also give its proof, using results which will 
be established later on. 

Theorem 3.1. Let tp £ ^T 1 (1R 3 ) ; with \\tp\\ = 1, and ^jv(x) = IljLi t P{ x j)- Denote by ipN,t the 
time evolution of ipN determined by the Schrddinger equation jl-4\) > an d let J^ tl for k = 1, . . . , N, 
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be the marginal distributions corresponding to ipN,t- Assume that the coupling constant A > —4/tt. 
Then for all fixed iGl, we have, for N — > oo ; 



fc>i 



r N ,t = {7$}*>i - r t = { 7i (fc) } fc >i on c 1 = 4 

respect to the product of the weak* topologies of C\. Here 

k 

where (ft is the solution of the non-linear Hartree equation 



(3.1) 



(3.2) 



id m = (1 - A) 1 / 2 ^ + A (J- * |^ t rt 



(3.3) 



Urti/i initial data (ft=o = <~p- 
Remarks. 



i) The factorization assumption for the initial state V'iv is only necessary in the limit N — > oo. 
In other words, if 7^ denotes the fc-particle marginal associated with tp^, it is enough to 
assume that, for every fixed k > 1, 



Tr 



as AT — > 00 



where 7^ c )(x/ c ; x' fc ) = 1X7=1 V 9 ( x j)v( x j) f° r some 95 G -ff^lR 3 ). 

ii) The result can also be easily generalized to the case where the initial state is not in a pure 
state. If we assume the initial density matrix to be given by 

N 

7at(x;x / ) = nV 1 ^;^-)' 
3=1 

for some 7M S -ff^M 3 ) x H l {R 3 ), we can prove the convergence (j3.1|) ; in this case the limit 

r t = {7* }fc>l is given by 7 ^(x fe ;xJ t ) = 11^=1 7* ( x ji x 'j) where 7^ (sc; ac') is the solution of 
the nonlinear equation 



i9t 7 (1) 



(l-A) 1 /2 + A (i*p t ) 



with pt(x) = 7 ( - 1 * ) (x; x) 



with initial data 7^ = 7W. This equation is equivalent to the Hartree equation ()3.3|) . if 

7 ( ^ 1 )(x;x') = Lp(x)Tp(x'). 

hi) Although we state our theorem specifically for the Coulomb potential V(x) = A|x| _1 , it is 
clear that our result and our proof still hold true for V(x) = A(x)\x\~ 1 +B(x), if, for example 
A,B € S(M. 3 ), the Schwarz class. 
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Proof. To simplify the notation we will assume t > 0. The proof for t < is similar. 

Step 1. Introduction of the cutoffs. For k > we define the regularized version of the initial 
one-particle wave function (p € i/ 1 (M 3 ) by 

ip K = exp(-K,\p\/N)(p 

with p = — iV. Note that the regularized wave function (p K depends on N. We use the notation 

N 

for the N body wave function corresponding to the regularized (p K . We denote by ip^t its time 
evolution with respect to the physical Hamiltonian Hjy (jl.lj) . and by i^jif t its time evolution with 

respect to the modified Hamiltonian defined in (|1.8|) . We also use the notation J^ Kt and 

7^ = 7^^. £ t for the family of marginal distributions corresponding to tp 1 ^ t and ipf} e t , respectively. 

Step 2. A priori bounds. In Sectional Theorem 14.41 we prove that for any k > 1 and for all 
e, k > there exists Nq = No(k, e, k) such that 

Tr \S 1 ...S k j%ls k ...S 1 \ <C k (3.4) 
for all N > Nq and for all t € M. The constant C is independent of k, e, k, N, t. 

Step 3. Compactness ofT^,t- Fix T > 0. From Theorem 15.11 it follows that, for every fixed 
k > 1, the sequence 7^ G C([0, T], 7^.) is compact with respect to the metric p k (defined in 

(|2,3|0 . This implies that the sequence = {^\} k=1 is compact in the space C([0, T],H) = 
®k>i T], TL k ) with respect to the product of the topologies generated by the metrics p k . Here 
we need the following standard argument (choice of the diagonal subsequence) to prove that, given 
a sequence Nj, with Nj — > 00 for j — > 00, there exists a subsequence Mj of Nj with Mj — > 00 
for j — > 00, such that 7^ t converges as j — » 00, for all fc > 1; by the compactness of 7^ there 

is a subsequence aij of Nj, with aij — * 00 for j — ► 00, such that 7^ . t converges. Next, by 

~(2) 

the compactness of r f Nt , there exists a subsequence «2j of «ij with Q?2.jf — * 00 and such that 

7alj t converges as j — > 00 (clearly, since «2,j is a subsequence of aij, 7^ t converges as well). 
Inductively we can define subsequences ol^h for all ell > 1, and we can set Mj = ctjj: then Mj — > 00 
for j — > 00 and 7^ t converges as j — > 00, for all > 1. 

It also follows from Theorem 15.11 that, if = {7<^t}fc>i £ C([0,T],H) denotes an arbitrary 
limit point of IV t, then we have 

Tr \S 1 ...S k -y£] t S k ...S 1 \ <C k (3.5) 
for all t € [0, T] and for all k > 1 (with the same constant C as in (|3.4[> ). 

Step 4- Convergence to the infinite hierarchy. In Theorem 16.11 we prove that any limit point 
Too,* £ C([0,T],Ti.) of the sequence T^t (with respect to the product of the topologies p k ) is a 
solution of the infinite hierarchy (|1.6[) . with initial value r oo ^ = o = {7o }fc>ij where 

k 

^ k \x k ;x' k ) = y»(a; i )<p(a^) (3.6) 
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(the K-dependence of ip K disappears in the limit ./V — > oo). 

Step 5. Uniqueness of the infinite hierarchy. In Theorem 17. 1 1 we demonstrate that for any given 
To = {To }fc>i £ Hi satisfying the estimate 

\h { k) \\ Hk =Tr \S 1 ...S k ^ k) S k ...S 1 \ <C k 

(k) 

for some C > 0, there is at most one solution = {7oo 4 }fc>i € C([0,T],H) of the infinite 

hierarchy (|1,6|) which satisfies r oo,t=o = To an d 

h£% k < C k (3.7) 

for all t € [0, T] and k > 1. This, together with (|3,5|) . and with the fact that, by Step 4, every 
limit point of T jv,t satisfies the infinite hierarchy (jl.fi j) , immediately implies that the sequence T jv,t 
converges for iV — > oo (a compact sequence with only one limit point is always convergent). Next, 
we note that the family of density T t = {7$ }ft>i, with 7^ defined by ()3.2j) . satisfies ()1.6j) and the 
initial condition ()3.6j) (the existence and uniqueness of a global solution of the nonlinear Hartree 
equation (|3.3j) . for A > A cr it, is proven in |2j). Moreover, it also satisfies the bound (j3.7|) . In fact, 
by Lemma 19.31 

|| 7 ? ) \\n k = Tr |Sj . . . S ki Ps k ... 5x| = fo, (1 - A) 1 /^)* < C fc . 

It follows that 

rjv,t -> T 4 as N ->• 00 

on C([0,T],H), with respect to the product of the topologies induced by p k , for fc > 1. This, in 
particular, implies that, for every fixed t G [0,T], we have rjv.t — ► T t in with respect to the 
product of the weak* topologies of TL k (which, since the sequence 7^ is uniformly bounded in N, 
are equivalent to the topologies generated by the metrics p k ). In turn this implies that, for every 
fixed t E [0,T], we have 

f Ntt ^Tt on ^ = 04, (3.8) 

k>l 

with respect to the product of the weak* topologies on C\, for k > 1. 

.Step Removal of the cutoffs. Recall that ipN,t denotes the evolution of the initial function 
■i/Tv generated by the Hamiltonian Hjy, defined in (jl.lj) . Moreover, ip% t and ipff t denote the time 
evolution of the regularized initial wave function ip 1 ^ with respect to the dynamics generated by the 
original Hamiltonian Hn and, respectively, by the regularized Hamiltonian Hjy (defined in (|1,8|0 . 
We remind that 7^, l^\ t and 7^ = T^ K et denote the marginal distributions associated with 
ipN,t, to 4>Nt anci > respectively, with tp^f t . In Proposition l<S. 1| we show that 

||^-V^II<c^ 1/4 

for every N large enough (depending on e, n). This clearly implies that 

Tr|7^ K -7{f,Ul<^,V4 (3 . 9) 

for all k > 1 and for all N large enough (depending on e, k). Moreover, in Proposition 18.21 we show 
that 

W N ,t -1>N,t\\ = Un-^nW <Ck 
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for all N and t G R. Therefore 



for all k > 1, t G R and N > 1. 



Tr|7^L-7Sl<^ 



(3.10) 



Step 7. Conclusion of the proof. For every fixed k > 1 and £ G [0, T], and for every compact 
operator G JC k: we have 



Tr 



(fc) (fc) 
7jv,t " 7t 



< !|J (fe) ||TV 



(fc) _ (jfc) 
7jv,t 7iv,K,t 



+ ||J (fc) ||Tr 



(k) _ (k) 



+ 



Tr 



7t 



where 7^ is defined in (|3.2j) . By (|3.9() and (|3.1U() . for any fixed S > 0, we can find k and e small 
enough such that both the first and the second term on the r.h.s. of the last equation are smaller 
than 5/3, uniformly in N, if N is large enough. Finally, by (|3.8|) . the last term can be made 
smaller than 5/3 by choosing N sufficiently large (and keeping e and k fixed). Thus, for every fixed 
t G [0,T], k > 1, G K k , and 5 > 0, we have 



Tr J« f 7 « 



(fc) 

7oo,t 



< 5 



for all N large enough. In other words, Tjvt — ► with respect to the product of the weak* 
topologies of C\, k > 1, for every fixed t G [0, T\. Since T < oo is arbitrary, this completes the 
proof of the theorem. □ 



4 A— Priori Estimates 

The aim of this section is to prove an priori bound, of the form 

Tr \S 1 ...S k jW t S k ...S l \ <C k , (4.1) 

for the marginals 7^ associated with the solution tp^ K t of the modified Schrodinger equation (|1.11|) 

with regularized initial data V'jv^o = ( recan that Sj = (1 — A.,-) 1 / 4 ). To prove (|4.1j) we derive 

an upper and a lower bound for powers of Hn in terms of corresponding powers of the kinetic 
energy. 

Proposition 4.1 (Lower Bound for H^). We distinguish two cases. 

i) If A > 0: For every k > 0, and C < 1, there exists Nq = No(k,C,e) (depending also on the 
cutoff e in the potential) such that 

(V, E%1>) > C k {N k (4>, Sf... S^) + N h -\^ SfSi . . . S 2 k _^)} (4.2) 

for all N > Nq, and for all ip invariant with respect to permutations. 

ii) If — 4/-7T < A < 0: For every k > 0, and C < 1, there exists Nq = No(k, C, e) such that 

H%^) > C k (l + -A)^ {N k (^, Sf... S&) + N k ~\^ SfSi . . . S 2 k _^)} (4.3) 

for all N > Nq, and for all ij) symmetric with respect to permutations. Here [k/2] = k/2 if k 
is even, and [k/2] = (k + l)/2 if k is odd. 
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Proof. We use a two-step induction over k. For k = the statement is trivial. For k = 1, and 
A > 0, the claim follows by the positivity of the potential. For k = 1 and A < 0, it follows from the 
operator inequality 

1 <"S 1 2 



|xi — X2\ 2 

(see Lemma 19. II part (i)) and by the symmetry of if) with respect to permutations. Next we assume 
that the claim holds true for k = n, and we prove it for k = n + 2. We consider the case A < (if 
A > the proof is easier), and we assume that n is even (for odd n the proof is analogous). By the 
induction assumption, for any C < 1, there exists No(n,C,e) such that 

ty, K +2 ^) = (V, H N H n N H N ^) > C n N n (l + ^A)"/ 2 ty, H N S 2 . . . S 2 n H N ^) (4.4) 

for all N > No(n, C, e), and for all ip symmetric with respect to permutations (note that also H^ip 
is symmetric). Here we are neglecting, because of its positivity, the contribution arising from the 
second term in the parenthesis on the r.h.s. of (|4.3[) . Writing 

N n \ N I 

h n = j2 s] + h N with h N = J2 S ] 



N ^— ' \xi — xA + eN' 1 

j=n+l j=l Kj 



equation (j4.4j) implies that 



^H^)>C n N n (l + ^Af/ 2 (ip,S 2 Sl..S 2 n S 2 iP) 

i,j>n+l 



+ C^N n (l + jXT /2 £ {(?P, S 2 S 2 ... S 2 n h N iP) + c.c.) , 



IT . 

i>n+l 

where c.c. denotes the complex conjugate. Using the permutation symmetry we find 



(tP,H% +2 4>) > C n N n (l + ^A) n/2 



(N — n)(N — n — Sf... S^) + (N - S{S' 2 . . . S^) 

+ ^-f;^ ^, VI.SI • • • S 2 n+ ^) + c.c. (4.5) 
+ A (^-^y — l) ( ^ ^ . . . ^ + c . c . 
(Ar- n )(Ar_ n _i)(Ar_ n _ 2 ) 2 2 

+ ^ — w, v n+2tn+3 s 1 • • • s n+1 ip) + c.c. 

where we use the notation Vfj = l/(\xi — Xj \ + eN' 1 ). We first estimate the last term (for positive 
A, this term is positive, and thus can be neglected). Since V^ +2 n+3 commutes with Sj, for all 
j < n + 1, we have 

K+2,n+3^i 2 • • • Sn+iVO = (ip,S 1 ... S n+1 V^ n+3 S n+1 . . . S^) < | (V, 5? . . . S 2 +2 V>) (4.6) 

because V^ +2 n+3 < \x n +2~ x n +3\ _1 < ( 7r /2)»S' 2 + 2 fLemma l9.ll part (i)). Now we turn our attention 
to the fourth term in (|4.5|) . We have 

|(^^n+25l 2 ...S 2 +1 ^)| = |(V,5 n+1 ...5 2 ^ n+2 S 2 S 2 ...5 n+1 ^| 

= |O0j 5'n+2'S'n+i • • • S±{S 1 S n+2 Vi tn+2 SiS n+2 }Si . . . S n+ iS n+2 V>)| (4.7) 
< const (V>, Sf ... S 2 +2 %ji) 
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where we used that, by Lemma l9~Tl part (ii), the norm \\S t S n+2 ^in+2^^n+2 II * s finite, uniformly 
in e. Finally we consider the third term in (|4.5|) . To this end we note that 

(if>, V{ )2 Sf . . . S 2 n+ ^) = (v, s n+1 . . . s 3 vi 2 sfsls 3 . . . S n+1 V0 

= (V>, S n+1 . . . S 3 S 2 S 2 1 {S^ 2 S 2 1 Vl 2 S 2 }S 2 1 S 2 . . . S n+1 i;) . 

For every 5 > there exists, by Lemma |9~T1 part (hi), a constant C £i $ > (depending also on e) 
such that 

| (V, Vl 2 Sl . . . S 2 n+l i,) I < C £:6 N s (iP, S\S 2 2 . . . S 2 n+l ^) . (4.8) 
We fix < 5 < 1. Then, inserting (|4~o% (|4~T)) . and (jHHJ) into (j!3|) . we find 

> + JA)^ {(iV - n)(iV - n - 1) + ^A (iV = n){N - U ~ ^ N - U - - 2) 

- const — K^,5f... S n+2 ip) 

+ C "iV"(l + ^A)"/ 2 | (JV — n) — C £i5 7V 5 } ty, SfS 2 . . . S 2 n+1 ^) 
= C"N^ 2 (1 + - x) n/2(N-n)(N-n-l) 



x (l + ^A^^ - const 2±i) (^, 5? . . . S 2 +2 ^) 



Next, since C < 1, we can find No(n + 2, C, e) > Nq{u, C, e) such that the four inequalities 

(N-n)[N-n- 1) _ /„ % . N — n — 2 n + 1\ /„ vr \ 

> C, 1 + -A const — — > C ( 1 + -A ) and 



- ' V 4 N N ) ~ \ 4 

*-» 2C , ri-c,^i> c 



iV ~ ' V N-n. 
are satisfied for all A r > A r o(n + 2, C, e). This implies that 

(V, i^+ 2 vO > c-+ 2 (i + 1 a) (n+2)/2 {N^, S 2 ... S 2 n+2 ^) + iV« +1 (^, SfS 2 . . . S 2 n ^)} 

for all N > No(n + 2, C), and completes the proof of the proposition. □ 

In the next proposition, we prove an upper bound for powers of H^. 

Proposition 4.2 (Upper Bound for H^). For I > 1, we use the notation ate = (aj, . . . ,ag) £ 
N+, and |oe^| = oji + ■ • • + ote- Assume C is large enough (depending only on X). Then, for every 
fixed k > 0, e > and for every a > there exists < C(k, e, a) < oo (also depending on X) such 
that 

fc-l 

{ij,H k N ^)<C k ^,{S 2 + ...S 2 N ) k ^) + C{k,e,a)N a Y J £ N'^S**...^) (4.9) 

€=1 CX£:|<X{ |=fc 

/or a// N and for all ip symmetric with respect to permutations of the N particles. 
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Proof. We prove the proposition by a two step induction over k. The statement is clear for k = 
and k = 1 (in this case one can choose C(k = l,e, a) = 0). We assume now that the statement 
holds for k = n, and we prove it for k = n + 2. From the induction assumption, with a replaced by 
a{n — l)/n, and since tp and H^ip are symmetric with respect to permutations, we have 

(i,,H n N +2 ^) = ^,H N H n N H N ^) 

< C n (?p, H N (S 2 + ■■■ + S 2 N ) n H N ^) 

( 4 -!0) 

+ C{n,e,a{n-l)/n)N< n ~ 1 ^ n Y J N r (iP,H N S 2ai . . . S? ar H N iP) . 

r=l oi r :\a. r \=k 

We start by considering the first term on the r.h.s. of the last equation. Applying a Schwarz 
inequality and using again the permutation symmetry, we find 

(^,H N (S 2 + --- + S 2 N ) n H N ^) 

< 2(V, (S 2 + --- + S 2 N ) n+2 ^) + N(N - 1)A 2 (V, Vl 2 {S 2 + ■■■ + S 2 N ) n Vl 2 ^) . 

The weighted Minkowski inequality 

(S? + • • • + S 2 N ) n < 2(5| + • • • + Sjf) n + c(n){Sl n + S 2 2 n ) 

(where the constant c(n) can be chosen as c n n n , for some c > 0) implies that 

(V, H N (S 2 + ■■■ + S 2 N ) n H N iP) < 2 (V, (S 2 + ■ ■ ■ + S 2 N ) n+2 i;) 

+ 2N(N - 1)A 2 (V>, V{ 2 {S 2 + ■■■ + S 2 N ) n V{^) (4.11) 
+ AN 2 X 2 c(n)^,Vl 2 S 2n Vl 2 ^). 

The second term on the r.h.s. of the last equation can be bounded by 

2iV(iV - 1)A 2 (^, Vl 2 (S 2 + ■■■ + S 2 N ) n V{^) <CN(N- (S 2 + ... S 2 N ) n S 2 S^) 

< C (^(S 2 + --- + S 2 N r +2 iP) 



(4.12) 



because Vf 2 commutes with (S 2 + • • • + Sfj), and (V± 2 ) 2 < \x\ — x 2 \~ 2 < C S\S 2 (see Lemma l9~Tjl 
and because N(N — 1)(V>, S 2 S 2 ip) < {ip, (S 2 + • • • + Sjj) 2 ip) for permutation invariant ip. Next we 
consider the last term on the r.h.s. of ()4.11|) . Using Lemma 19.21 for every a > 0, n > and e > 0, 
we find a constant D(n, e, a) such that 



iV 2 Vl 2 S\ n Vl 2 < D{n, e, a)N a N 2 £ sf n - m) SfslN^ 1 . 



-,2{n—m) c 4 c>2 ; 
m=l 

In the sum over m, we can add more factors of S 2 , so that the order of each monomial is 2{n + 2) 
(and we can compare this contribution with the sum in the r.h.s. of (|4.9|) ). More precisely we can 
bound 

q2(n-m) q4 q2 _ q2(n-m+2) Q 2 ^ Q2(n-m+2) ^2 q2 q2 q2 
°l >3 1I D 2 — D 1 <3 2 ^ &i D 2 D 3 . . . D m D m+1 

because S 2 > 1. Hence, for every a > 0, there exists a constant C\{n, e, a) (also depending on A) 
such that 



n+l 



4c(n)X 2 N 2 (^Vl 2 S 2n Vl 2 ^) < d(n, e, a)N a £ £ ^ ty, S 2 " 1 . . . Sf ^) 

£=1 ct^:|c*^|=n+2 
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From the last equation, together with (|4.11[) . ()4.12|) . it follows that 
(V, H N (Sf + ■■■ + S 2 N ) n H N ^) < C 2 (i;, (S? + • • • + S 2 N ) n+2 ij) 



n+l 



+ C 1 (n, £ ,a)iV^ {^S 2 ^...S 2 rm l ^ 

l=\ OLf.\oLf\=n+2 

provided C is large enough (depending only on A). 

Next we consider the second term on the r.h.s. of (|4.1Uj) . Applying the Schwarz inequality, we 
find 

n-l 

Y, E N " W> 5n S " ai ■ ■ S 2 r ar H N ^) 
r=l a r :\oL r \=n 

n-l 

^ 2 E E N r (i;,(S 2 + --- + S 2 N ) 2 S 2a \..S 2a ^) (4.14) 

r=l a r :\a r \=n 
n-l 

+E E Yi^vfA ai ■ ■ ■ s r ar n^) ■ 

r=l a. r :\oL r \=n i<j 

Using the permutation symmetry, the first term on the r.h.s. of last equation can be rewritten as 
n-l 

E E N " W> (sf + --- + S 2 N ) 2 S 2 ^ . . . S 2a ^) 

r=l ot r :\ot r \=n 
n-l 

= E E N "' { r S l (ai+2) Si"' ■ ■ ■ + 2r ( r " X ) Sf ai+1) Sl [a2+l) Sl a * . . . S 2a ^) 

r=l a r :\a r \=n 

+ 2r(N - r)ty, S? ai+1) . . . S 2 ^S 2 +1 ^) + (N - r)(if>, S 2 ^ . . . S 2 ^S^) 
+2(N - r )(N-r- S 2ai . . . S 2 ^ S 2 r+1 S 2 +2 ^)} 

(4.15) 

which clearly implies that 

n— 1 n+l 

E E Nr w> (sf + --- + s 2 N ) 2 s 2 ^ . . . s 2a ^) < e E Nt w, sf ai ■ ■ ■ sf v) • 

r=l a r :\a r \=n £=1 <y.g : \ ct( | =n+2 

The second term on the r.h.s. of (|4.14|) . on the other hand, can be rewritten as 

n-l 

E E ^ r E(^^A 2ai ---^^) 

r=l a r :\a. r \=n i<j 

E E ^{^^^(^,^^...5^^) (416) 



r=l a r :\a. r \=n 

+ r(N - r)ty, Vl r+1 S 2 r • . . S^V^) 
+- 



.1 (at - r)(N - r - Vf +1>r+2 S 2a > . . . S 2a ^ +W i,) 
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To bound the last contribution we use that, since V^ +l r+2 commutes with S\,...,S r , and since 
(ip, V r e + i^ r+2 S 2ai . . . S 2oir V^ +ltr+2 tp) < C (ip, S 2ai . . . S 2oir S 2 +l S 2 +2 %l)) ■ 



v r+l,r+2 — ° J r+l J r+2) 



(4.17) 



Hence 

n-l 

E E ^(«-r)(JV-r-l)(«W + 2?...? r W 

r=l OL r :\oL r \=n 

n+l 

^ C E E n'ms^.-.s^). 

i=\ a.^:\cy.(\=n+2 

As for the second term in (|4.16|) . we write 

(V, ^r+i ^ 2ai • • • S^Vl r+l i>) = ty, S*» . . . ^ r+1 S 2ai ^ r+1 ^ . . . S 2 Q » 
and then we observe that, by Lemma 19.21 for every a > we can find D(n, e, a) such that 

ai 

Vtr+i S i aiV hr+i < D(n, e, a)N a l n E 5 2{oi " m+2) 5 r 2 +1 N^ 1 . (4.18) 

m=l 

This implies that, for a suitable constant C2(n, e, a), 
n-l 



E E ^M-^ - r){4>, Vf tr+1 Sl ai . . . Sr ar Vi jr+1 ip) 

r=l ol t : | OLr | =n 

n— 1 ai 

< C 2 (n, e, a) iV a /™ E E E ^ r+m (V>, S 2K - m+2) S 2 2a2 . . . S 2 .^ S 2 r+1 ^) (4.19) 



n+l 

<C 2 (n,e,a)Ar a / n E E iY '^' S ^ ' ' ' 5 ?» ■ 

i=l ai :\at£ |=n+2 

In the last inequality we used that in the sum over m we can add more derivatives (until the total 
order of each summand is 2(n + 2)), according to the bound 

iirf+m p2(ai-m+2) n2a 2 n2a r q2 ^ Ajr+m Q 2(ai-m+2) Q 2a 2 o2o r q2 c<2 q2 

iV O-^ 2 • • • ^v-\-\ — 1 2 ' ' ' T 7*4-1 r+2 ' ' ' r-\-m' 

Moreover, note that r + m < r + at\ = r + n— (a 2 + • • • + at r ) — r + n — (r — 1) < n+l (this 
explains why the sum over the index £, in the last line of Q4.19J1 . only goes up to n + 1). 

Finally, we consider the first term in (|4.16f> . We have, using the symmetry with respect to 
permutations, 

(V, V£ 2 . . . sf^v^) = (V, S£ 3 . . . s^vi 2 s 2 ^s 2 2 a wi 2 s^ . . . S^) 

< (v, s%*... s^vi 2 s 2{ai+a2) vi 2 s^ . . . 5°». 

Applying again Lemma 19.21 we find, for every a > 0, a constant D(n,e,a) such that 
Vl 2 sf ai+a2) Vl 2 <D{n,e,a)N a l n E S 2{ai+Q2 ~ m+2) S^N^ 1 . 

m=l 
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Hence 

n—l 



E E ^(r-l)^,^^ 1 ...^^) 

r=l QL r :|tx r |=n 

n—l ai+«2 

<C 3 (n,e,a)iV»/"E E E N^-\il>,^ ai ^-^ 2) £%S** j>) (4.20) 

r=l ct r :[ct r |=n m=l 
n+1 

< C 3 (n, e, a)N a l n E E AT^, S 2 ^ . . . Sf ^) . 

l=\ tx£-.\a£ |=n+2 

Similarly to (|4.19|) . in order to control the terms in the sum over m, we used the trivial bound 

Arr+m-l q2(ai+a2-m+2) Q 2 q2o 3 Q 2a r ^ AT r+m-l Q2(oi+02-m+2) Q 2 c2a 3 r-2o r q2 q2 

iV *_)-^ *^2 3 * ' ' v — 1 2 3 ' * * v r+1 ' * * v-\-rri 1 ' 

and the inequality r + m — 1 < r — 1 + a\ + at2 < r + k — 1 — («3 + . . . a r ) < n + 1. 
Inserting (|CT7jl . (jCTIj) . and (j4~2llj) in (j1~TK)) we find 

n-i 

E E ^ r E^ y M- 5 i 2ai --- 5 ' ar ^) 

r=l o: r :|ct r |=ra «<j 

n+1 

< C 4 (n, £ ,a)iV a /"E E JVV,Sf 1 ...SfV). 

£=1 cx£:|<X£ |=n+2 

Combining the last bound and (|4.15f) with Q4.14JI we find 

n-l 

E E N r ^,H N Sl a K..S 2 ^H N ^) 

r—l ct r :\oL r \=n 

n+1 

< C 5 (n, e, a)N a / n E E N ' ^ ■ ■ ■ " 

£=1 a^:|a^|=n+2 

Defining C(n + 2, e, a) = C n C\ (n, e, a) + C(n, e, a(n — l)/n) • C5 (n, e, a) , the last equation, together 
with (|4"?T3|) and (JUTUJ), implies that 

< C n+2 (V, (5? + • • • + S%T^) 

n+1 

+ C(n + 2, a, e)iV a E E ^ ^ ^ ' ' • 

£=1 c^:|c*^|=n+2 

and completes the proof of the proposition. □ 

In order to apply the last proposition to our initial state Vat we need the following lemma, 
which shows, together with Proposition 14.21 that 

(r N ,H k N r N )<c k (4.21) 

for N large enough (depending on e, k and k). 
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Lemma 4.3. Fix ip € ^(R 3 ) with \\ip\\ = 1. For k > de/me ip K = exp(-K\p\/N)ip. Let 

^(x) = nf=i^(^)- 

^ VKe /iaue ||v3 K || < 1 and 

fa*, (1 + p 2 ) fc /V) < C fc (l + fc! (N/K) k - 2 )\\<pf Hl . 

ii) For every k > 1 and k > i/tere exists a constant C(k, n) such that 

k-i 

J2 E N e (r N ,Sf a K..S 2 e a ^ N )<C(k,K)N k - 1 . 

£=1 cxe:\a e \=k 

Hi) For fixed k, and k, there exists No(k,n) such that 

(r N ,(s 2 + ...s 2 N ) k r N )<c k N k 

for all N > Nq (the constant C is independent of k, k and N). 
Proof, i) The inequality \\f K \\ < 1 follows from e~ K \ p \/ N < 1. Next we compute 

(<p K , (1 + p 2 ) k / 2 <p K ) < 2 k l 2 J dp (1 + \p\ k )e- 2K ^ ,N \(p{p)\ 2 

<2 fc / 2 ^l + sup^pl^e- 2 ^!/^^ J dp(l + \p\ 2 )\<p(p)\ 2 

<C k (l + k\(N/ K ) k - 2 )M 2 Hl . 

ii) For any £ < k — 1 and ct£ = (ai, . . . , at) G N+, with |a^| = k, we have, using part i), 

e 

N'wh, s 2ai . . . s 2a *r N ) < n* H(<p k , (i + p 2 r^ 2 ^) 

3=1 

I 

< N k N e ' k Y[ C Qj (l + ctj\ K 2 -^N^- 2 )\\ip\\ 2 H1 

3=1 

e 

< C k N k \\<p\\ 2 £ YliN 1 -^ + aj ! K 2 ~ a i N- 1 ) . 

3=1 

Since ctj > 1 for all j = 1, . . . ,£, and since £ < k — 1 there is at least one j <G {1, . . . ,£} such that 
«j > 2. Thus, 

fc-i fc 

£ E W'tyfr, ^i 2Ql • • • Sf>ft) < C*IMIni W*" 1 JJ(1 + fe!K 2 " fc ) fc = N k ~ 1 C(k, k) . 

i=l OLf.\a t \=k 3=1 

iii) By the permutation symmetry of ip^, and by part ii), we have 

fc-i 

(^,(5 2 + ..- + ^) fc ^)<iV fc (^^i 2 ...^) + E E A^,S^...S 2 ^) 

<=1 c^:|a^|=fc 

< N k {<p K , (1 + p 2 ) 1 / 2 ^)^ + C{k, k)^- 1 

< C fc iV fc f 1 + ^^iv- 1 
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where the constant C only depends on For fixed k, k we can now choose iVo large enough 

such that C(k, k)N^ < 1. Then for N > N we have (ifj%, (S? + • • • + S 2 N ) k ^) < (2C) k N k , which 
proves the claim. □ 

Using the bound (|4.21|) . the conservation of the energy, and the lower bound for iJ^ given in 
Proposition 14. 1| we finally arrive at the a-priori bound l|4.1jl . 

Theorem 4.4 (A-Priori Estimate). Fix a one-particle wave function (p € ff 1 (IR 3 ), with ||(£>|| = 1. 
For k > let tp K = exp(— K\p\/N)ip. Define ^>jy-(x) = Y\f=i f K ( x j)- Let ip 1 ^. be the solution of the 
Schrddinger equation hi. 11)) with modified interaction and initial data ip 1 ^- Then, for all fixed 
e, k > and k > 1 there exists Nq(e, k, k) such that 

(^J, Sf . . . Sl < C k and (^£ t , SfSi . . . Sg_i < C k N (4.22) 

~(k) 

for all N > Nq and for all t G K. The constant C only depends on A and \\ip\\ui. Denoting by y Nt 
the k-particle marginal distribution associated to V'yvt' i4-^>$ * s equivalent to 



Tr 



Si...St y£j St - ..Si 
for all N > N . 

Proof. From Proposition 14.21 we have 



< C k 



and 



Tr 



Sf #2 • • • Sk-i y^\Sk-i ■ ■ ■ S2S 2 



k-1 



< C k N 



(^ N ,H k N r N )<CHrN,(Sl + --- + S 2 N ) k rN) + C(k,e,a)N a Y; E NttyhS?* 1 . 

By Lemma l4,3| part ii) and iii) we find 

(Vfr, H%^) < C k N k + C(k, e, k, a)N a N k - 1 < C k N k (l + c ( k '^ a ) N a-i 

for N large enough, depending on k and k. Hence, fixing a < 1, we find 

U%,HUN)<CfN k 

for all N large enough, depending on k, e, re, and with a constant C\ only depending on and 
A. By the conservation of the energy and by Proposition 14.11 we get 

c k N k > (r N , H k N r N ) = w%% H k N ^£ t ) 

> C k N k {^ v Sl... Sl + C k N k -\^ v SfS 2 . . . SU fift) 



for all N large enough. 



□ 



5 Compactness of the sequence Vnj. 



Recall that we defined Fjv,t = {^N Vk^i as ^ e marginal densities associated to the solution ip^ K t 
of the regularized Schodinger equation with initial data ip 1 ^. They satisfy the regularized 

BBGKY hierarchy 

k 



:1 = E[( 1 - A ^) 1/2 ^S]+ A ^ 1 E [v°te-*i) 



(*•■) 



~{k) 

TNt 



3=1 



l<i<j<k 



(5.1) 



+ A ( 1 "] v f)E Tr fc+i V%x 3 -x k+l )^ t 



5=1 
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where V £ (x) = (\x\ + eN~ 1 )~ 1 . 

Theorem 5.1. Fix k > 1 and T > 0. Then 7^' . G C([0,T],7-^) /or a// Af /arge enough (depending 

on k, e and k). Moreover the sequence 7^ is compact in C([0, T],TL k ) with respect to the metric 

p k (defined in and, if 7^ G C([0, T],7ifc) denotes an arbitrary limit point of the sequence 

(fc) 



1 00, t 

7]^ j wra£/i respect to the metric p k , we have 



h£] t \\H k = Tr\S x ... S kl t] t S k . . . 5x| < C k (5.2) 
for all t G [0, T] (^/ie constant C is independent of k,e, and k). 

~(k) . . ^ 

Proof. We prove that the sequence "f Nt is equicontinuous in t with respect to the metric p k . To 
this end we rewrite the BBGKY hierarchy (|5.1I) in the integral form 

k „+ , k 



~(fc) ~{k) 
TN t = TN s ~ 1 



3=1 Js i<3 Js 

p. . 

1 ** 8 



\ ~(fe+i)i 



where Trfc + i denotes the partial trace over the k + l'th variable. Next we choose G K. k such 
that ^SjJ^'Sj || is bounded for all j = 1, . . . k, and we compute 



Tr( fc ) J« (7^-7^) = I ^ {k) {STWSj - S^Sri) S^S, 

j= i Js 

"4^ /^^^ {sj 1 J {k) S 3 Sj l V £ {x l -x,)S^S^ T S, 
i<j Js 

-S J jWSf 1 Sf$j T S j Sf 1 V'(x i - x^Sj 1 ) 

~ <(! " aT^E fdrTr^J^S^Vixj - x k+1 )S£ v S k+1 j^ T S k+1 ] 
3=1 Js 

where we used the fact that S k+ i = (1 + p 2 k+ i) 1 ^ 4 commutes with j( k \ Since we assumed 
IjS'j J^Sy 1 1| to be bounded for all j = 1, ... k and since \\SJ V £ (xi — x^S^W < C (see Lemma 
19. we find, using Theorem 14.41 



7k t - v 



<c fc |t 



for a constant which depends on k and JW, but is independent of t,s G [0, T] and of N, for 
iV large enough (depending on k, e and k). Since the set of all j( fc ) G /C fc , with ||5 j jWS'^ 1 || 

bounded for j = 1, . . . k, is a dense subset of A k , it follows by Lemma I9~4l that the sequence j N t 
is equicontinuous with respect to the metric p k (defined in (|2.2j0 . This in particular implies that 
7^j G C([0, T],H k ) for all A" large enough. Moreover, since the sequence 7^ G C([0,T],Ti k ) is 

uniformly bounded, it follows by the Arzela-Ascoli Theorem (see jllj . Theorem 1.28) that it is 

(k) 

compact with respect to the metric p k , and that any limit point 7^ satisfies (|5.2|) . □ 
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6 Convergence to the infinite hierarchy 

In this section we consider limit points of the sequence Tnj = {7jvt}fcLi> where are the 
marginals associated with the solution ip^ K t of the modified Schrodinger equation 1)1. lip with reg- 
ularized initial data ijj 1 ^. We prove that any limit point of this sequence, as N —* oo, is a solution 
of the infinite hierarchy of equations (|1.6p . 

Theorem 6.1. Suppose = {7<2\}fc>i £ C([0,T],7i,) is a limit point of the sequence Tj^ >t = 

{7^l}fcLi with respect to the product of the topologies induced by the metrics p^ (defined in \2. C J\) ). 
Then T^ t satisfies the infinite hierarchy hl.b)) with initial data 

k 

7So(x fc ; x' fc ) = 7( f } (x fc ; x' fc ) = ip(xj)<p(xj) . (6.1) 

i=i 

Proof. Without loss of generality we can assume that r./v,t — > Too,* G C([0,T],7i) with respect to 
the product of the topologies induced by the metrics pk- Then, for every fixed k > 1 and t G [0, T], 
we have 7^ — > 7^ € Wfe with respect to the weak* topology of Hk (because the sequence 7^ 
is uniformly bounded in TLk and the metric p^ is equivalent to the weak* topology for uniformly 
bounded sequences). Hence, for every t £ [0,T], k > 1, and every S .Afc, we have 

Tr J (fc) (7$ " 72!*) - foriV^ TO . ( 6 . 2) 

We choose j( fc ) 6 /Cfc such that SjJ^S^ 1 is bounded for all j = 1, ... fc, and we obtain, from the 

__4 rT j 

regularized BBGKY hierarchy (|5.ip . rewritten in integral form, 
Tr (fc) = Xr (fe) j(fc) W (fc) (t) ^) o 

-iA(l r^Tr^ 1 ) J«wM(t-*) [^ fe+1 ,7^ 1) ] 
i= i ^0 

where 

W«(i) 7 (fc) = e- itH{k) ^e UH{k) with F( fc ) = E 5 1 + ^E V i,3 

j=l i<j 

and where we use the notation Vfj = V £ (x{ — Xj). Next we consider the limit N —* 00 of ()6.3p . 

From (|6.2p . the l.h.s. of ()6.3p converges to Tr <7^7qJ>V ^ s ^ or ^ ne ^ rs ^ term on the r.h.s. of (|6.3I) 
we find 

TrV J« (^)(t)7f -W «(t)7So) =TVWJ« (w<«(t)-^(t))^ 

+ Tr« J«^(t) (7^0 -7Sq" 

where we used the notation 



(6.3) 



(6.4) 
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for the free evolution of the first k particles. The first contribution on the r.h.s. of (|6.4|) can be 
handled as follows: 



and thus, using the permutation symmetry of , and the a-priori estimate from Theorem 14.41 

TrW j(fc) ( w (*)( t ) 7^1 < CiV-^^HJ^H ||^ 2 S' 1 - 1 S' 2 - 1 ||Tr( fe ) SiSs^^i 

< CA" 1 

with a constant C depending on k, on the observable and on i. As for the second term on the 
r.h.s. of Q6.4J1 we have 

TrWjWwf(t)(7« -7S0) =Tr« (l<«(-f)J«) (7$ - 7 £o) "> 

for A -> 00, because, if G /C fc , then also U^\-t)jW G £ fc C A, and thus Q can be used. 
Next we prove that the difference 



l-k/N)J2 fdsTr^ jWllW(t-s) [V^,^] 



j= l JO 



converges to zero, for A — > 00. Here and henceforth, Vi = V(xi Xj J — I X 2 X j .To this end, 
we rewrite it as the sum of four terms 

k . 

-k/N^ I d^ {k+1) ^W W (t)[^ +1 ,7g] 
3=1 J o 

3=1 JO V J 

k . 

+ £ [ dsTr^ jMU?\t-s)[Vf !k+1 -V jik+ in^ 

■ 1 J 



3=1 J0 



+ £ / d«Tr fc+1 J^U$\t - s) [ Vj , k+ i, (7$, - 7 &>. 
j=l J 

The first term converges to zero, for A — > 00, because 

fc/A^ Td^Tr^) ^^(t)^,^] 



< CF 1 ^ Wl^S^S^W sup Tr( fc+1 ) |5i5fc+i 7 g 

se[o,t] 

< C^N- 1 sup Tr( fc+1 ) S 1 S k+1 -y^ 1) S k+1 S 1 < C^N" 1 

se[o,t] 
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for a constant C k t depending on k, on t, and on the observable J^ k \ To control the second term 
we note that 



o 



i=i 



Ak+l) 



k k 



JO 



t-s 



-iXN^ 1 E E ( ds /„ d ^ Tr(fc+1) J {k) €\t -s-r) \V e * m M k Hr) \V[ k+1 ,t. 

j=l £<m 



(k+1) 



Writing down the four terms arising from the two commutators and using the permutation sym- 
metry, we get the bound 



E I dsT * (k+1) J(k) (uW(t-s)-4 k \t-s)) [Vf^tf* 
< C^k^N- 1 sup 1 11 c r(k)a-l 



(k+1) 

s 



i<k 



+ 



x sup (Tr |5i5 2 7^ 1) 5 3 | + Tr |SiS 2 S 3 7]£ 
se[o,T] 



(k+i) 

s 



< C^kN- 1 / 2 sup Tr S 1 S 2 Stf%+ 1) S 1 S 2 S 3 <CN- 1 / 2 
se[o,T] 

where the constant C depends on e,k,t and on j( k \ Here we used that HS^^Zmll < Ce~~ l / 2 N l l 2 , 
the symmetry of 7^^ 1 \ and Theorem 14.41 To bound the third term in ©, we note that 



V j,k+1 ~ V j,k+1 



eN- 1 



\xj - x k+1 \(\xj - x k+1 \ + eN- 



-IV 



Hence 



E / dsTr™ J^U^(t - s) \v£ k+1 - V jik+1 ,^ 
3=1 







< kteN- 1 H^fV^^ill l^f 1 



1 



jrS x 1 S 2 2 1| sup Tr \SiS 2 'y'fi^'~' S: 



(k+1) , 



\xi - x 2 \(\xi - x 2 \ + eN x ) 



se[o,T] 



< Cfc^iV- 1 / 2 sup Tr S 1 S 2 ^ 1) S 1 S 2 < ^ _1/2 



E jf d,sTr fc+1 jWl4 k \t-s) \v jM1 , (^f -7^ 1} 



(6.5) 



se[o,T] 

for a constant depending on k, t, and J^. Finally we consider the last term in @ 
k <-t 

This term converges to zero, because, since £ /Cfc C Afe, we have (Uq (s — t)J^)Vj tk +i € .4fc+i 
for every j < k, every s and In fact 

iisr 1 • • • s k iM h) (* - 1)^)^!^ 1 • • • sSiii 
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This proves that the integrand in Q6.5|) converges to zero as N — ► oo, for every s G [0, t] and for 
every j = 1, . . . , k. Since the integrand is uniformly bounded in s, it follows that (|6.5[) converges 
to zero for N — ► oo. 

We have proven that, for every fixed t G [0, T], k > 1 and G Kk such that SjJ^Sj 1 is 
finite for all j, we have 



Tr« jW 7 g t = Tr J< fc >W< fc >(thS> " Td.Tr^ 1 ) J«Z<«(t _ s) \ Vj , k+h ^] . (6.6) 

'o 



Since the set of G K,k such that sup J<fc < oo is a dense subset of .4.^, it follows by 

a simple approximation argument, that ()6.6|) holds true for all G Ak- Thus 



for all t G [0, T]. Finally we have to prove that 7^ is given by (|6,1() . Recall that 



7Si(x fe ;xl,) = n^(^)^(4) 

with ip K = exp(— n\p\/N)ip. Hence, for any G ICk, we have 

TV J« ( 7 f - 7So) = Tr J« ( 7 « - 7^0) + Tr J<*> ($> - 7oo ,o 
The second term converges to zero, for N — » 00. As for the first one, we have 



(*) 



(fc) ~(fc) 



< Ch\\(p- <p K \\ < CkN- 1 



(6.7) 



TrJ« ( 7 « -7$) I <IU (fc) II Tr 
where we used that ||v? K || < \\<p\\ = 1, and that 

h-<p K f = J dp(l-e- K ^ N ) 2 \ip(p)\ 2 <*>N- 2 \\<p\\ H i. 

Since the choice of -/V on the right side of Q6.7JI is arbitrary, we find 

TVJ« (7f-7So)=0 

for all j(® G /Cfc. Moreover, using the fact that 7 q and 7^ nave finite Tif. norm, a simple 
approximation argument shows that the last equation is true for all G Ak- This proves that 



(fc) (fe) 
7oo o = 7o • 



□ 



7 Uniqueness of the solution of the infinite hierarchy 

In this section we prove the uniqueness of the solution of the infinite hierarchy Ijl.fij) . We already 



know that Tt = {7^ }&>!, with 



7f ) (x fc ;x' fc ) = [J<pt(a;j)¥>t(a 
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and tpt the solution of the nonlinear Hartree equation (jl.3|) . is a solution of (jl.6|) . Since, by Theorem 
16.11 we know that every limit point of T^,t is a solution of (|1.6|) . it follows that Tt is the only limit 
point of rjv,t- 

Theorem 7.1. Fix T > 0. Let T = {7 fe) } fc >i G H = fe >i?4- Moreover, suppose that there 
exists C > suc/j i/iaf 

||7? } ||H fc = Tr |5i . . . srf^fc • • • Si | < C k 

for all k > 1. TTien i/iere exists at most one solution T t £ C([0, T],H) o/ i/ie infinite hierarchy 
11. b)) such that T t =o = and 

hl k) \\H k <c» 

for all t € [0,T]. 

Proof. Suppose Tt = {7t^}fc>i is a solution of the infinite BBGKY hierarchy Ijl.fijl . so that Tt=o = 
Tq and 



(*)i 



\H, 



(7.1) 



for all t G [0, T] and > 1. Rewriting (|l.(i|) in integral form we find 



It 



1 -uS*\t)^-i\J2 [ t dslt§'\t-B)Tr» 1 \y jM u'ti k+ % (7-2) 

7 = 1 J ° 



where we use the notation Vij = V(xi | 1 and where the free evolution UQ k \t) is 

defined by 

W (fc) (t) 7 (&) = exp ^-itJ^S^j 7 W exp ^*X>1 
We iterate (|7.2|) . and find 



7 P =ui k \t) 1 < £ ) + ^ C(M,"0 + C(M,n) 



ra-1 



m=l 



where 



(7.3) 



fc fc+m— 1 

Cftt,m) = (-ar^... £ / d Sl 

.(fc+m-1), 



ds m U^ k \t - st) Tv k+ i Vj lik +i,Uo'^ XJ ( s i ~ s 2 



/(*+!)/ 







(*m-l S m ) Tr/j^ m Vj mj fc_|_ m , Z^q (^m)To 



(fe+m)/ \ (fc+m) 



and the error term n, t) is given by 

k k+n— 1 „j 

C(M,n) = (-iA)*J>- £ / ds i- 

il = l jn = l J ° 



ds n ^ fe) (i - si) Tr fc+ i Vj-j jfe+ i , U]f +1> (si - s 2 



/(fc+i), 



x Z^j (s n _i - s n ) Tr fc+n Vf n ,fc+n>7s 



(fc+n) 



(7.4) 
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In order to bound the error term, we note that, for any i > 1, j = 1, . . . £, s € R, and 7^ +1 ) e T~Lt+i, 
we have 



\\u^\s)^ i+ i[v jMl ^]\\ Hl 

= Tr® \Si... S t U®{a) (Tt t+1 [V j/+1 , 7 (m) ]) S i ...S 1 \ 
< TrW \Sx . . . S e (tt £+1 V j4+ i 7 {m) ) S t . . . S x \ 
+ TrW |5i . . . S< (Tr m 7 ^ +1 ) K- j+i) S e . . . 5i| 



(7.5) 



because the free evolution ZYg commutes with the operators Sj. We consider the first term. By 
the cyclicity of the partial trace and because 

Tr^\Tr i+1 A\ <Tr<- t+r >\A\, 

(see Proposition 9.4 in 0]), we have 

TrW \SL-.Se (Tr m V jt m7 (m) ) 5/. . . 5 X | 

= TrW [5! . . . S e (Tr e+1 S^V^j^ S i+1 ) S t . . . Si| 

< Tr^ +1 ) |5i . . . SeS^Vu+nV+VSt+tft . . . Si| 

= Tr^ +1 )| (SjS^Vj^S^S^ Si... S i+l7 ^S i+1 ...Si\ 

< WS^V^S-^W || 7 (m) lk +1 < C|| 7 (m) lk +1 

for a constant C, independent of In the last inequality we used Lemma 19. 11 part (ii). The second 
term on the r.h.s. of ()7.5j) can be bounded similarly. We get 

\\U^(s)Tr l+ i[V j>e+ i,^]\\ Ht < 2C\h^\\ He+1 ■ 

Applying this bound iteratively to (|7.4|) we find 

\\C(k,n,t)\\n k <2 n X n C n t n ( k + n ) sup || 7 ( fe+ri ) 

V n J ae[0,t] 

By 1)7. 1|) we have, 



"Hk+n - 



\\((k,n,t)\\<D k D n t n , (7.6) 

where D depends only on the constant C in the bound (|7.1|) (and on A). Now suppose that 
Fit = {7it }fc>i an d T2t = {72t }fc>i are two solutions of the infinite BBGKY hierarchy with 
ri^ = o = r2,t=o = To and satisfying l|7.1|) . Then, for f < 1/(2D), we have 

hi? -7£ } |k = llCi(M,n) -b(k,t,n)\\<H h < 2D k 2- n 

for any n > 1 (note that the terms C(fc, i, m) in the sum over m in (|7.3[) depend only on the initial 
data 7g fe ^ , and thus do not contribute to the difference 7^ — 7^ ) . Since n > 1 is arbitrary, we find 

(fc) (k) 
Ti,t = 72,t 

for all k > 1 and for all t < 1/(2D). Since the bound (|7.1|) holds uniformly in i, for t € [0, T], the 
argument can be iterated to prove that Fij = F2,t, for all t G [0, T] . □ 
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8 Removal of the cutoffs 

From Theorem 16,11 and Theorem 17. 11 we know that r^t converges to F t = {7t*^}jfc>i G C([0, T\,TL) 
for iV — > oo, with respect to the product of the topologies pk (defined in 12, 3 j) . In this section, 
we show how to remove the cutoffs e and k, which are used to regularize the interaction and the 
initial data. To remove the cutoff e we need to compare two dynamics, the one generated by the 
modified Hamiltonian l|1.8|) . and the one generated by the original Hamiltonian (jl.lj) : this task is 
accomplished in the next proposition. 

Proposition 8.1. Suppose ip^t ^ s a solution of the modified Schrodinger equation 11.11)) with 
initial data i/j^-, an d let V'jvt be the solution of the original Schrodinger equation \l.<j\) > with the 
same initial data. Then there exists a constant C, independent of N , k, e and t such that 

u K N e ,-r N ,t\\<cte^ 

for all t > 0, and for all N large enough (depending on e and k). 
Proof. We use the notation <^v,t = V'jvt an d 4>N,t = "^jvt- We compute 

dt \\<l>N,t - (pNjW 2 = i{HN<t>N,t ~ H N <fr N j, 4>N,t ~ <t>N,t) ~ i(<t>N,t ~ 4>N,t, H N (f>N,t - H N <fr N j) 

= -Im (0jv,t - <j>N,t, (Hn ~ H N )(j)N,t) 
by the self-adjointess of H]y. Hence 

±d t Un,t ~ M\ 2 < \\</>N,t ~ j>N,t\\ \\(Hn ~ H N )0 N , t \\ . (8.1) 

We have 

H _ 5 = Af (_L 1 A _ a f sn-i 

N^\\xi-Xj\ Ixi-xA+eN- 1 ) N ^ \xi - xA(\xi - xA + eN- 1 ) 
Therefore, using the symmetry with respect to permutations 
\\{H N - H N )0 Ntt f = ($N,t, (H N - H N ) 2 N , t ) 

< X 2 E 2 (<j)N,t, 1 FT] i— m f" 77TTT </>JV,i) 

Fi - a^KFi - x 2 \ + eN L ) \x s - X4\(\x3 — x±\ + eN L ) 

+ A 2 e 2 iV _1 $ Nit , , | m , $N,t) 

\xi - X2\(\xi - X2\+ eN 1 ) |x 2 - 2C3KF2 - £3 1 + J 

+ AViV^ (^, ki _ X2|2(|xi l X2|+gjV _ 1)2 fet) • 

(8.2) 

The first term on the r.h.s. of the last equation can be estimated by 

1 1 



' \xi - x 2 |(|xi - x 2 | + eN x ) \x 3 - x 4 |(|x 3 - x 4 | + eN 

<c$ N , t , sls 2 2 slsl4> N>t ) < c 



< {^N,ti 1 FTi fn 4>N,t) 

Fi - F3 - ^4F 
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for all t G K, and for every N large enough (depending on e, k, see Theorem 14. 4(1 . The second term 
can be bounded, using again Theorem 14.41 by 

(0N,t, 1 17| j— Ar _iN 1 m r i rTTTT «AlV,t) 

Fi - »2|(fi - %2\ +eN 1 ) \x 2 — xs\{\X2 — x s \ + eiV i J 



< e 1 N(^ N ,t, i rrn r 0Ar,t) 



(8.3) 



Fl ~~ x 2| 2 F2 - S3 1 

< Ce- 1 iV(^,53 1 L_^ 3 ^ t ) 

Fi — x 2| 

^e- 1 ^, S 2 S 2 S 2 ^ N , t ) 

< Ce~ l N 

for all N large enough. Finally, we estimate the last term on the r.h.s. of (|8.2|) as follows 

N-^l^-^l+.iV-^ ^ ^ ^ /2iV3/2 (^ 1^-^5/2 ^) 

< Ce-WN 3 ' 2 $ NtU Sf /a sffc) 
by Lemma part (i). Next we note that, by Theorem 14.41 

&N,t,SfS$$ N , t ) < $N,t,StSi0 N>t ) < CN 

for all N large enough, and for all t (the constant C is independent of e). Since (S1S2) 1 / 2 < 
N 1 / 2 + N- 1 / 2 SiS 2 , we find 

@ N ,t,sl /2 $ /2 M < NV 2 $ N>t ,sM<t>N, t ) + N-v 2 $ Njt ,s!si$ Njt ) < cn 1 / 2 . 

From (jH.Hj) we get 

for all N large enough. From (|8.2|) it follows that 

\\(H N -H N )$ N £ <Ce 1 / 2 

for all N large enough (depending on e), and for all t (the constant C only depends on A and 
||</?||i/i, see Theorem 14. 4[) . Hence, by (|8.1(l and by Gronwall's Lemma, 

±WN,t -0N,tf< Ce l/i U N ,t ~ fell => life ~ fell < Cel ' H 

for all t > and for all N large enough. □ 

Finally we have to remove the cutoff k from the initial wave function tjj 1 ^. 

Proposition 8.2. Suppose (p G i? 1 (M 3 ), with \\tp\\ = 1. For k > put ip K = exp(— K,\p\/N)ip. 
Suppose V'jvt anc ^ V'Af.i are the solutions of the Schrodinger equation \1.J\) with initial data ip^(x) = 

Art-fell < K\\<p\\ H i 



Ylj=i f K { x 3); an d, respectively, ^jv(x) = 11^=1 r/ien 



/or a/Z t G R and all N > 1. 
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Proof. By the unitarity of the time evolution, we have 



^N,t-^N,t\\ = I^AT-^VII 



N 



N 



3=1 3=1 

N (8.4) 
< Yl H^ K ( x i) • • • ^fai-iX^fe) - (p(xj))<p(x j+ i) . . . <p(x N )\\ 

3=1 

<N\\<p«-<p\\. 

Here we used that the L 2 norm of ip K is bounded by \\<p K \\ < \\<p\\ = 1. Since p K = e~ K ^ p ^ N (p and 
1 _ e -K\p\/N < kN' 1 \p\, we have 



dp {e-^ /N - l) 5 ''\(p{p)\ 2 < k 2 N- 2 j &p\p\ 2 \(p{p)\ 2 = k 2 N 



:-' _ ,.2 a — 2 i _l|2^ i 



By (|Q|) . we find 



□ 



9 Some Technical Results 

In this section we collect technical results used throughout the paper. 

In the first lemma we show how to control singularities like \x\ — X2| _< \ for a < 3 in terms of 
the operators Sj = (1 + p 2 ) 1/4 = (1 - Aj) 1 / 4 , for j = 1, 2. 

Lemma 9.1. Let V(x) = | a; | 1 and V £ (x) = (\x\ + eiV~ 1 )~ 1 . Moreover, we set Vij = V(xi — Xj) 
and Vfj = V s (xi — Xj). 

i) For all a < 3 there exists C(a) such that 



1 



- < C(a) S?S% for all a, f3 > with a + 13 = 2a. 



\xi - x 2 

If a = 1, we have the tight bound 



Vi >2 < ~ S?S$ for all a,(3>0 with a + (3 = 2 . 



ii) The operator 

is bounded. Moreover 



S\S 2 1 Vi^S 1 1 S 2 1 



uniformly in e. 
Hi) For all 5 > 0, there exists C$ < oo such that 

WS^S^Vf^W < C s e- 5 N 5 



(9.1) 
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Proof. To prove i) we use that 

^—<^Sf, and | l —<C(a)S? a (9.2) 

\X\ — X2\ 2, \X\ — X2\ 

for some constant C{a) < oo, and for i = 1,2. This is proven in [H]. The statement i) can now be 
shown by the following general result. Suppose A > 0, and B, C > 1 are three operators, with 

i4 < 5 and A < C 

and such that B commutes with C. Then 

A < B a C fS for all a, (3 with a + (3 = 1 . (9.3) 

In fact, A < B implies that B~ l l 2 AB~ l l 2 < 1, and, by the operator monotonicity of powers smaller 
than one, this implies that {B~ x l 2 AB~ l l 2 ) a < 1 for all a < 1. On the other hand A < C implies 
B- l / 2 AB- 1 / 2 < B- l C and also (B~V 2 AB- 1 / 2 )? < B'^C^ for all < 1. Hence 

3-1/2^-1/2 = ^ B -1I2 AB -XI2^I2^ B -XI% AB -1I^^ B --XI-1 AB -XI2^I% 

< (B-^AB-^f < B^C 13 . 

Multiplying both sides by B 1 ! 2 gives flOJ). Part i) follows by (fP)) and by taking A = 

\x\ — X2\~ a , B = S\ and C = S2. Next we prove ii). We will make use of following fact (known as 
the Holmgren-Shur inequality). Suppose that the set of mutually orthognal projections {P n }^ =1 
resolves the identity in a strong sense, that is 



N 

A — v 



lim V P n = 1 . 



n=l 

Then for any operator A we have the norm inequality 



\\A\\ <max sup V ||P fc AP n || , sup V ||-Pfc^-Pn|| • (9.4) 

V k n=l n fc=l / 

We choose 

P fc := ^(le^ 1 < v\ + K 16") , 
then the boundedness of Si^^Vi^S*! -1 ^ 1 follows by (|9.4|) if we prove that 

ll^fe Si 1 S 2 ' -, r S 1 S 2 ' P n \\ <C 2-l fc -«l/ 2 . (9.5) 

\X\ — X2\ 

for a constant C independent of n, k. To show the estimate (|9.5|) we consider two possibilities: 
k > n and n > k. 

When k > n we can bound the l.h.s. of (|9.5[) by 

HP^r 1 !! U^ 1 1 r S^|| H^PJI . (9.6) 

\Xi-X2\ 



Since, by part i 

II ^9 

\X\ - X2 



Sry ~ " Sex II 



2 U. „.| "2 II — 2 ' 
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and 



\PkS^\\ < 2- k+ \ ||5iP n || < 2 r 



we obtain (|9,5|) (in fact a stronger estimate). If n > k we write 
1 



Pk Si S 2 



\X\ — X 2 



Si S 2 P n 
= P k Si 1 S 2 1 ( P 2 1 + l)(pi + l)- 1 



1 



|xi - x 2 



Si S 2 Pn 



PkS^S^ipl + l) 



\xi - x 2 



+ p fc 5r i s 2 - i (pf + i) 



(pi + -0 , 



i 



|»i - x 2 | 



5 1 5 2 - i P n . (9.7) 



The norm of the first contribution can be now bounded by 

\\P k St\\\\S 2 - 1 r -!—-S 2 - 1 \\\\Si 3 P n \\ < C2- 3 l fc -"l, 

\Xi X 2 \ 

hence we are left with the task of checking that also the norm of second contribution in (|9.7|) 
satisfies a similar bound. Note now that 



(pi + 1)- 1 



i 



\xi - x 2 

therefore it suffices to estimate 



(pi + 1)- 1 pi 



X\ - x 2 



+ 



Xl - X 2 



\xi - x 2 \ 3 \Xl - x 2 \ 

p i- s ~ ' - s 7 ' /" ' rr—, % (pi + i)' 1 s i s 2 l PJ 



pi (p\ + 1)- 1 



\Xi - X 2 \ 



and 



\P k Si L S 



1 c-l xi - x 2 



2 \x\ - a?2| 3 
We bound (|9.8|) . using the result of part i), by 



2 M IIC-3 , -1.,,, . X l 

\x\ - x 2 \ 3 



\\P k S(\\\\Si"S 2 ' Pl 
As for ()9.9j) we estimate it by 



P^pi + iy 1 SiS^Pni 



S^S^WWSi 2 P n \\ < C2- 2 \ k ~ n \. 



(9.8) 
(9.9) 



\ Pk s{' 2 \\ \\sf 2 s ? . pi s-^ 5r 5/2 n ii5r 1/2 p n 



Fi - ^2 



1 



1 



1151 S2 \x 1 - X2 \^ uu \x 1 -x 2 \m 



Si 1/2 s 2 l \ 



< C2~ lk - nl/2 



hence the result. The same line of reasoning can clearly be applied also to prove (|9.1|) . uniformly 
in e > 0. Finally, we show part iii). To this end, we write 



/•oo o— 2 o 

Si 2 S 2 Wl 2 S 2 = Si 2 Vl 2 + J&S s-^ * (p 2 • + V^ 2 • g + i 2 + 
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The first term is bounded, by part i), uniformly in e. As for the integral, for any given 5 > its 
norm can be bounded by 



(e- l N) 5 



ds 



o s 



1/4 



\p2 



s + 1 + p l 2 



S 



1 



1 \xi-x 2 \ 



1 



1 



+ 



(lp 2 l + l) 1/2 - 



8+1 + p\ 



\x\ -x 2 \ x ~ h (|pa| + l) 1 - 6 
1 



s 2 (N + i) 



l-<5 



s + 1 + p A 2 



(ipai+ir^sr 3 - 



^1 — X2| 3 / 2 



IX1-X2IV2 (|p2| + 1)V2 

r°° 1 



S2\P2\(\P2\ + 1) 1/2 



S + l +p\ 



< C s e- 5 N 5 . 



□ 



The following lemma is used in the proof of Proposition 14,21 in order to control terms of the 
form V^ 2 Si n Vi 2 by powers of S\ and S 2 . 

Lemma 9.2. For every n > 1, e > 0, and < a < 1, there exists D(n,e,a), independent of N, 
such that 



2(ra-m+2) g2jym-l 



m=l 



Vl 2 S\ n Vl 2 < D(n,e,a)N a £ S 2 
Proof. Suppose first that n = 2i is even. Then 

S 2n = Sf = {l+p\f<2 i {l + (pi)") < 2 e + 6 e J2 Pi 



2( 



a=l 



where pi >a denotes the a's component of the vector p\. Hence 

Vl 2 Sf V{ 2 < 2" {Vl 2 f + 6^ Vl 2 p\* a V{ >2 . 



(9.10) 



a=l 



Using the commutator expansion 



Vl 2 p{ a = p{ a V£ 2 + Ptj [• • • [Vl&Pl,a],Pl, a ], ■ ■ -Pl,a] 

i=i W * „ ' 



j commutators 



and a Schwarz inequality, we obtain the bound 



Vl 2 < 2p{ a (Vl 2 fpl a + 2*J2 I [■■■M, 2 ,Pi,a], ..-Pi 

3=1 ' : • / -~ — 



1 2 1-3 
I P\A 



(9.11) 



j commutators 



Next we use (V^) 2 <CS 2 S% and 

I [••• [Vl 2 ,Pl,a],---Pl,c 



< C(j\) 2 



j commutators 



< C(j\)\e- L N) 



(|xi -x 2 \ + eN~ 1 ) 2 i+ 2 



1 Nr\ 2 3~ 1 + a . 



\Xl - x 2 



3— a 



<D 1 (j,e,a)N 2 i- 1+a SfS 2 
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for some < a < 1 (the constant Di(j,e,a) is proportional to (j!) 2 , to e 2j+1 a , and diverges 
logarithmically in a, for a — > 0). From (|9.11l) we find 

Vl 2P i%Vl 2 < C P { a SMp{ a + D 2 (£,e,a) j^pf^ SfS 2 N 2 ^ 

for a new constant D 2 (£,e,a) independent of N. Inserting the last equation into (|9.1U|) . summing 
over a, and using again (V^ 2 ) 2 < CS 2 S 2 , it follows that 

Vi^SfV^ < D(£, e, a) ^Sf i+2 S 2 + £ s} (t - j) St% iV 2 ^ 1+a j . 
Replacing the indices £ and j by re = 2£ and m = 2j, we have 

^ 2 5i 2n ^ 2 < C(n,e,a) £ S 2{n ~ m+2) S 2 2 N m ~ 1+a 

m=l 

which proves the claim if n is even (since all terms in the sum over m are positive, we can also 
allow m to be odd). 

If n is odd, the proof is a little bit more difficult. Let n = 2£ + 1. Then 



^ 2 sf +2 ^ 2 = ^ 2 (i +piys(v 1 \ 



2\£ c2t/£ 



< 2^ 2 S?V?; 2 + 6< £ F^S 2 ^ 

< C{£) U(^ 2 ) 2 ^ + |[^ 2 , 5 X ]| 2 + Y.PlM V l2) 2 SiP{« +pU[Vi% S 1 }\V 1) , 

V a=l 

3 

+ E E I [• • • W;a.Pi,«], • • • iPi,a] \ 2 s 1P {J 



a=l j=l 

3 e 



j commutators 



+ E E <a IK; • • [^2,Pl,a] • • ^iiiViri 

a=l j=l ^ 



j commutators 



(9.12) 



Next we note that, for any integer r > 0, 

3 



9 r 



• J poo 

E/ d«. 



r+l 



s + l+p 2 \ dxi^dx\ 



-Vf 2 + h.c. 



1 



s + l+pf 



2 ■ 



Using 



Qr+lye 



< c 



(r + l)! 



(|asi -x 2 \ +eN- l Y+ 2 ' 



dxi^dx\ ( 

we find, for every r > 0, e > and a > 0, a constant C(r, e, a), such that 

0' 



dx r u 



— 2 q— 1 
1 °2 



< C(r,e,a)N r+ ^ 
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and thus, for every a > 0, we have the operator inequality 

2 



d r 



< C{r,e,a)SfS%N 2r+a 



This can be used in the second, fourth and last term on the r.h.s. of 1)9. 12j) . The other terms can 
be handled as in the case of even n. □ 

In the next lemma we give a proof of the fact that a solution ipt of the Hartree equation 1)3.3}) has 
i? 1 / 2 -norm uniformly bounded in time. We assume here that the initial data p G H 1 / 2 (we apply 
this result in the proof of Theorem l3.11 Step 5, where we have the stronger condition p € H 1 ( 



Lemma 9.3. Suppose p G /f 1//2 (IR 3 ), with \\p\\ = 1, and let pt be the solution of the nonlinear 
Hartree equation 

id t p t = (1 - A)^ 2 p t + A(|-r * \pt\ 2 )pt 

with initial data pt=o = P- Then, if A > — 4/-7T, there exists a constant C , depending only on A 
such that 

{p t ,{i-^) l i 2 p t )<c\\p\\ m/ , 

for all t € R. 

Proof. The L 2 norm of p is conserved, so \\pt\\ = 1- Also the Hartree energy 

E(p)= [ dx\(l-A) 1 /Mx)\ 2 + ^ [ dxd Vr ^— \p{x)\ 2 \p{y)\ 2 
J 2 J \x-y\ 



is conserved by the time evolution. Note that 

1 



dxdy • 



\x - y\ 



\f( x )\ \f(y)\ < su p / d y 



< 



7T||V|| 



r^l^)l 2 lbll 2 

\x - y\ 

dy\(l-A y )^p( y )\ 2 



where we used the operator inequality \x — y\ 1 < (7r/2)(l — A y ) 1//2 for every x € 
we find 

(1 + j\)(p, (1 - A) l ' 2 p) < E(p) < (p, (1 - A) 1 /^) . 
For A > 0, on the other hand, we have 



For A < 0, 



(p, (1 - A^p) < E{<p) < (1 + -X)(p, (1 - A) 1 ^) . 

Hence, for all A > — 4/V, we have 

(p t , (1 - A) l l 2 p t ) < CE(p t ) = CE{p) < C\\p\\ 2 m/2 
for a constant C only depending on A. 



□ 



Finally, in the next lemma, we give a criterium for the equicontinuity of a sequence of time 

(k) 

dependent density matrices j Nt E TCk, with respect to the metric p^- This result is used in the 

proof of Theorem 15 .1| in Step 3, to show the compactness of the sequence 7^ t G C([0, T], Hk) with 
respect to the metric p^- 
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(k) 

Lemma 9.4. A sequence of time- dependent density matrices j\r t , N = 1,2,..., defined for t € 
[0, T] and satisfying 

sup h$ t \\H k < C (9.13) 
te[o,T] 

/or a// iV, is equicontinuous in C([0,T],Tlk) with respect to the metric pk (defined in \2.ty) ). if and 
only if, for all in a dense subset of Ak, and for every rj > there exists a 5 > such that 



< rj 



(9.14) 



/or aZZ iV, whenever \t — s\ < 5. 



Proof. Equicontinuity w.r.t. the metric pk means that, for any rj > there exists 5 > (independent 
of N), such that 



, (fc) (fc) 



3=1 



rp r(fc) / (fc) (fc) 



< e 



(9.15) 



if \t - s\ < 5. Recall that {jf ] } j>i was chosen as a dense countable subset of the unit ball of Ak- 
Using (|9.13|) . one can approximate any E Ak by an appropriate finite linear combinations of 
the Jj k ^ and thus one can easily prove that 1)9.15(1 implies (|9.14j) . 

Next we prove the opposite implication. By a standard approximation argument, one can prove 
that, if flgJH) holds for a11 j(k) 

in a dense subset of Ak, then it holds for every J( fc ) in A k - In 
particular, for every j > 1, and n > 0, we can find 5(j, rj) > such that 



rp T (k) ( (k) (fc) \ 

Tr J) (7^ - 7^J 



< rj 



for all N, if |t — s\ < 5(j,r]). Moreover, using (|9.13|) . we note that, given r\ > 0, we have 



£ 2-i Tr Jf ( 7 W _ 7 W) | < £ 2-i|| Jf | Ufc (|| 7 W | k + || 7 « < C £ 2- < n/2 
if m is sufficiently large (independently of A?" and of i, s G [0, T]). Hence 



E 2 ~ 

i>l 



r(*) / ( fc ) ( fc ) 



< 77/2+^2^' 



T r( fe ) Z' ( fe ) ( fc ) 



With 5 = minj< m S(j, rj/2), we find 



< 



for all i,s £ [0, T] with |i — s| < <5 and for all N. This proves that (|9.14l) implies l|9.15|) . 



□ 
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